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Substantial Disparities Exist in 
Global Ocean Evaporation Estimates 

Global Mean Evap 60N/S (mm/d) 12 — mon smooth 







Observations" <-More Observ Constrained 


A Hierarchy of Global Ocean LHF Estimates 

AGCMs w/ Specified SSTs (AMIPs) GEOS-5, ERA-20CM Ensembles 

Incorporate best historical estimates of SST, sea ice, radiative forcing 
Atmospheric "weather noise" is inconsistent with specified SST so sfc fluxes can be 
wrong sign (e.g. Indian Ocean Monsoon, high latitude oceans). Averaging over 
ensemble members helps isolate SST-forced signal. 

Reduced Observational Reanalyses: NOAA 20CR V2C, ERA-20C 

Incorporate observed Sfc Press (20CR) and Marine Winds (ERA-20C) to recover much 
of the evaporation induced by true synoptic or weather "noise" 

Comprehensive Reanalyses (MERRA-2) 

Full suite of observational constraints- both conventional and remote sensing. 

But... substantial uncertanties owing to evolving satellite observing system. 
Multi-source Statistically Blended OAFlux, LargeYeager 

Blend reanalysis, satellite, and ocean buoy information. While climatological biases 
are removed, non-physical trends or variations in components remain. 

Satellite Retrievals GSSTF3, SeaFlux, HOAPS3... 

Global coverage. Retrieved near sfc wind speed, & humidity used with SST to drive 
accurate bulk aerodynamic flux estimates. Satellite inter-calibration, spacecraft 
pointing variations crucial. Short record (1987-present). 

In situ Measurements ICOADS, IVAD, Res Cruises 

VOS and buoys offer direct measurements. Sparse data coverage (esp south of 
30S. Changes in measurement techniques (e.g. shipboard anemometer height). 


Questions we'll address 


■Given the uncoupled framework of "AMIP" experiments, 
what can they tell us regarding evaporation variability? 

■ Do Reduced Observations Reanalyses (RedObs) using sfc 
pressure (and wind) provide a more realistic picture of 
evaporation variability? 

■What signals of interannual variability (e.g. ENSO) and 
decadal variability (IPO) are detectable with this hierarchy 
of evaporation estimates? 


SST variability forcing 
of evap is limited to key 
regions 


AMIP tightly couples 
evap to SST forcing 


Local Correlation (SST, Evap) 
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Correct extra-tropical 
SST forcing by evap is 
missing in AMIPs 


Strongest evap forcing 
of SST resides in 
extra-tropical eastern 


Local Correlation (SST f , Evap) 



ocean basins 


60N 

40N 

20N 

EQ 

20S 

40S 

60S 


I-20C 



60N 

40N 

20N 

EQ 

20S 

40S 

60S 


MERRA-2 


ERA-20CM 





_l I ZZd 

- 0.4 - 0.3 - 0.2 


- 0.1 


0.1 


I L 

0.2 0.3 



- 0.5 


0.4 


0.5 







What's the Interpretation? 


As with previous studies (e.g. Barsugli & Battisti,1999; Wu et 
al, 2006; Chen et al, 2013; 2014) specifying AGCM boundary 
forcing (i.e. SST) yields evaporation ("atmospheric noise") 
inconsistent with the SST forcing. 

- Atmospheric "weather noise" differs determ/n/st7cc///y from that which 
originally helped force the observed SSTs (even though it is 

stochasticcilly similar). 

- The problem is most pronounced at mid- to high latitudes and in 
association with the Monsoon regions. 

By ensemble averaging (Chen et al; 2013; 2014) we expect 
that evaporation associated with internal atmospheric 
variability is greatly reduced, leaving the SST forced signal. 


A Taylor Series Expansion of Bulk Aerodynamic 
Evaporation Around Monthly Climatology 


As in Richter and Xie (2008), Lorenz et al (2010) we first write the bulk 
formula for evaporation as a function of SST, Relative Humidity, Wind Speed, 


and Stability: 


E = C E p a Uq 0 [q s (SST)- RH • q s (SST + S) 



near-sfc q a 


where we use the analytical expression for saturation specific humidity 


q s (SST) = q/ (SST ] p = 

and q 0 , a, b, c are constants. 


S = T . -SST 

air 


Evaporation anomalies, SE, are expressed as 

8E = -^EdSST + ^ 8U + ^EdRH + ^ 8S + ^-<SC £ + res 
dSST dU dRH dS dC E E 

where he partial derivatives are "sensitivities" built from monthly resolved 
climatology and 5( ) denotes a monthly anomaly. 



Evaporation Change Mechanisms (AMIP Ensembles) 


Ocean 60 N/S time series (mm/d) 


12 mon smoothing 


GEOS5 AMIP Ensemble Mean 
Evap Decomposition Anomalies 
60N/S Oceans 12 — mo smth (r 
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ERA— 20CM AMIP Ensemble Mean 
Evap Decomposition Anomalies 
60N/S Oceans 12 — mo smth (mm/d) 



Rel Hum: 

Stability 

Residual 


1900 1910 1920 1930 1940 1950 1960 1970 1980 1990 2000 2010 


20 th Century Trends (% Climo Evap / Deg Global SST Change) 



ERA-20CM 1.17 5.59 -3.59 -0.50 -0.34 

•Actual evaporation trend lies substantially below C-C Rate due to 
offsetting contributions of increased RH, Stability and Exch Coeff. 
•Wind-related trends are small (but interannual signals are large). 








Evaporation Change Mechanisms (Reduced Obs Assimilation) 


Ocean 60 N/S time series (mm/d) 


12 mon smoothing 
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1900 1910 1920 1930 1940 1950 1960 1970 1980 1990 2000 2010 


20 th Century Trends (% Climo Evap / Deg Global SST Change) 



ERA-20C 4.95 5.50 -5.56 7.71 -2.70 


• Both NOAA 20CR and ERA-20C have more decadal variability than AMIPs 
•ERA-20C wind trend drives 20 th Century evaporation trend rivaling C-C rate. 








MERRA2 spec hum analysis increment z=1 (g/kg/d) 
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As SSMI contributes to near-sfc drying 
the RH decreases act to 
increase evaporation. 


1980 1985 1990 1995 2000 2005 2010 



1980 


F08 has major impact on 
GEOS5 wind speed bias 


2010 

Loss of SSMI suite 
but, scatterometer 
winds from ERS-2 
through 2011 and 
ASCAT on METOP-A 
continue to be 
assimilated. 


Reduced Observation Assimilation 

• Assimilating near-sfc weather (P s and/or wind speed) is a 
step toward recovery of evaporation (and other fluxes) 
driven by atmospheric state alone missed by AMIPs. 

• But... the inhomogeneity of data availability and its 
quality can induce significant uncertainties. 


Can we discern physical signals e.g. ENSO, PDO/IPO? 

ENSO: hi-pass (<7y) 1950/2010 data & composite by lags 
relative to Equatorial Pacific SST max 

IPO: Epoch difference (1999/2009 minus 1990/1999) 
fields. Period covers transition from warm to cold 
IPO phase (~ "hiatus") 


ERA-20C Composite 

El Nino Evaporation Anomalies (shaded, mm/d) 
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Lag (months) Relative to El Nino Peak 


EVAP Contributions: WSPD, QS-QA 
1.0a event composites 1950/2010 
Global Ocean (60 N/S) 


Does the Inter-decadal Pacific Oscillation 

Have a Robust Signal? 



Time series of IPO (Henley et al 2015) 


IPO Change 
Coincident with 
"Hiatus" 
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ERA-20C SST (°C) 1999/2009 minus 1990/1999 
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Inter-decadal Pacific Oscillation In ERA-20C 


• Consistent pattern of evaporation changes (shaded) relate to SST gradient 
patterns and altered low-level circulation. 

•Twin anticyclonic high pressure anomalies (green contours) develop in the 
eastern extra-tropical Pacific resulting from westward diabatic heating shift. 

• Off-equatorial wind speeds Eq to 20 N/S (black contours) consistent with 
surface pressure gradient changes enhance evaporation. 
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Summary Points: 


Consistent with coupled model results (e.g. Richter and Xie, 
2008; Lorenz et al, 2010), process related diagnostics indicate 
AMIP experiments show thermodynamic damping by RH and 
stability act to keep evaporation trends below C-C rate. 

Reduced Observations Reanalyses (RedObs) assimilating 
surface pressure (and wind) can provide a more consistent 
picture of evaporation variability. Some issues: wind speed 
shift @1940; dearth of So. Hem. observations. 

El Nino-related IA evap variations are large regionally. 
Equatorial wind speed decreases lead qs-qa maximum evap 
max tends to lag Eq SST max coherent global signals. 


Challenges: 


■ Data scrubbing / refurbishment a continuing need. Cross- 
comparisons among "hierarchy" products with varied data 
input are useful in ferreting out uncertainties. 

■Significant changes in data availability in the Satellite Era 
continue to challenge comprehensive reanalyses. More 
attention to reanalysis model analysis increments is needed 
to understand data / physics bias interplay. 

■ Future versions of satellite evaporation retrievals will likely 
have reduced trends due to better sensor calibration effects 
on wind speeds. Qs-qa biases need further study. Wx 
"regime" approaches a way forward. 
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